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We study the plasma-induced modifications of the potential and charge distribution across the 
interface of a plasma and a dielectric wall. For this purpose, the wall-bound surplus charge arising 
from the plasma is modelled as a quasi-stationary electron surface layer in thermal equilibrium 
with the wall. It satisfies Poisson's equation and minimizes the grand canonical potential of wall- 
thermalized excess electrons. Based on an effective model for a graded interface taking into account 
the image potential and the offset of the conduction band to the potential just outside the dielectric, 
we specifically calculate the modification of the potential and the distribution of the surplus electrons 
for MgO, Si02 and AI2O3 surfaces in contact with a helium discharge. Depending on the electron 
affinity of the surface, we find two vastly different behaviors. For negative electron affinity, electrons 
do not penetrate into the wall and a quasi-two-dimensional electron gas is formed in the image 
potential, while for positive electron affinity, electrons penetrate into the wall and a negative space 
charge layer develops in the interior of the dielectric. We also investigate how the non-neutral 
electron surface layer - which can be understood as the ultimate boundary of a bounded gas discharge 
- merges with the neutral bulk of the dielectric. 

PACS numbers: 52.40.Hf, 73.30.-|-y, 52.80.Tn 



I. INTRODUCTION 

O ^ Macroscopic objects in contact with an ionized g as ac- 
quire a negative charge because the influx of electrons 
from the plasma outruns the influx of ions. The col- 
^ . lection of electrons at the wall (boundary of the object) 
^»-^' gives rise to a repulsive Coulomb potential which reduces 
, the electron influx until the wall charge reaches a quasi- 
Q-f stationary value. As a consequence of the electron ac- 
] cumulation at the wall an electron depleted region, the 
plasma sheath, is formed adjacent to the wall. 

Most of the voltage driving the discharge drops across 
the sheath. Wall charges may however not only affect 
the spatial structure of the plasma but also surface- 
supported elementary process such as electron-ion re- 
combination and secondary electron emission, which 
are particularly important in dusty plasmasji"— dielec- 
tric barrier discharges;^"— and solid-state based micro- 
discharges j^"— A macroscopic description of the plasma- 
induced wall charge, sufficient for the modeling of the 
plasma sheath is clearly insufficient for quantifying the 
effect wall charges might have on these processes. A mi- 
croscopic description of the plasma-induced wall charge 
and the potential across the plasma- wall interface it leads 
to is required. 

Traditionally, plasma walls are treated as perfect 
absorbersJ^"— Irrespective of the microscopic interac- 
tion, all electrons and ions impinging on the wall are 
assumed to recombine instantly. From this model only 
the wall potential just outside the wall can be obtained. 
This is the potential that balances the electron and ion 
influx at the wall. A first, qualitative step going be- 
yond this model was taken by Emeleus and Coulter j^^d^ 
who envisaged the wall charge as a two-dimensional sur- 
face plasma coupled to the bulk plasma via electron- 
ion wall recombination. No attempt was however made 



to put this appealing idea onto a formal basis. Later 
the notion of a two-dimensional surface charge was de- 
veloped further by Behnke and coworkersii"— utilizing 
phenomenological rate equations for the electron and ion 
surface densities. In these equations, the microphysics 
at the wall is encapsulated in surface parameters, such 
as, electron and ion sticking coefficients, electron and 
ion desorption times, and an electron-ion wall recombi- 
nation coefficient. In principle these parameters can be 
calculated. Assuming, for instance, plasma electrons to 
ad- and desorb in the long-range image potential of the 
wall we calculated in our previous work electron sticking 
coefficients and electron desorption times for uncharged 
metallio^S and dielectric surfaces<^"— We also made a 
first attempt to estimate these two quantities for charged 
dielectric plasma walls^l and proposed a physisorption- 
inspired microscopic charging model for dust particles in 
a gas discharge 1^ 

In this work, we shift gears and focus on the potential 
and charge distribution across the plasma-wall interface 
after the quasi-stationary wall charge (electron adsorbate 
in the notation of our previous work^"— ) has been es- 
tablished. In other words, we extend the modeling of the 
plasma sheatbi^rii to the region inside the solid and cal- 
culate the plasma-induced modifications of the potential 
and charge distribution of the surface. Although knowing 
the potential and charge distribution across the interface 
may not be of particular importance for present day tech- 
nological plasmas, it is of fundamental interest from an 
interface physics point of view. In addition, considering 
the plasma wall as an integral part of the plasma sheath 
may become critical when the miniaturization of solid- 
state based plasma devices^"— continues. 

In the model outlined below we specifically consider 
a dielectric wall and treat the plasma-induced quasi- 
stationary wall charge, that is, the surplus charge on 
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top of the charge distribution of the bare, free-standing 
surface, as an electron surface layer of a certain extent, 
which is trapped by and in thermal equilibrium with the 
wall. In order to determine the chemical potential, width, 
and spatial position (relative to the crystallographic in- 
terface) of the electron surface layer, which depend on 
surface as well as plasma parameters, we employ a one- 
dimensional model for a graded interface between a col- 
lisionless plasma sheath and a dielectric surface which 
is assumed to be a perfect absorber, that is, the wall 
potential balances at a certain distance from the crystal- 
lographic interface the electron and the ion influx from 
the plasma. 

The model of a graded interface encompasses two im- 
portant ingredients: the surface dipole of the bare sur- 
face responsible for the offset of the conduction band 
minimum to the potential just outside the dielectric and 
the long-range image potential. The former accounts for 
the charge re-distribution of the free-standing, uncharged 
surface arising from the truncation of the crystal and 
the latter supports polarization-induced external surface 
states (image states), first predicted for liquid helium^^ 
and later studied for metallic and dielectric surface with 
negative electron affinity^Zr— , which may trap the elec- 
tron surface layer in front of the crystallographic inter- 
face. 

Originally proposed by Stern for the interface be- 
tween two dielectrics and later on used by others for 
semiconductor heterojunctions^ and electron trapping 
in nanoporeSf^ the graded interface model also guaran- 
tees continuity of the electrostatic potential across the 
plasma-wall interface. The model thus allows us to study 
the spatial distribution of the plasma-induced wall charge 
across the interface. To insert the surplus charge into the 
interface we follow Tkharev and Danilyuk^ and mini- 
mize, in the spirit of density functional theory,—"— the 
grand canonical potential of wall-thermalized excess elec- 
trons. We also investigate how the electron surface layer 
merges with the neutral bulk of the dielectric which we 
describe with the model of an intrinsic semiconductor. 

Various improvements of the model arc conceivable but 
the increased mathematical complexity would mask the 
general ideas we would like to convey. For instance, the 
model of a collisionless sheath could be replaced by more 
realistic modelsJ^"— Going beyond the perfect absorber 
model, on the other hand, is an unsolved problem. It 
would require the inclusion of electron desorption, elec- 
tron sticking, and electron-ion recombination, with the 
respective coefficients to be self-consistently calculated 
for the quasi-stationary electron adsorbate at the wall. 
Replacing the graded interface by an ab-initio theory 
for the surface, for instance, along the lines given in 
Refsj^i^i^i^, possibly taking ad-layers of the host gas' 
atoms or molecules as well as impurities and imperfec- 
tions into account, is desirable but at the present stage of 
the investigation impractical. It would require an expen- 
sive atomistic characterization of the plasma-wall inter- 
face, either experimentally via various surface diagnostics 




FIG. 1: Qualitative sketch of an interface between a plasma 
and a dielectric wall. Upper panel: Band structure, micro- 
scopic crystal potential merging with the image potential, 
and sheath potential. Lower panel: Effective potential for 
the graded interface on which the model of an electron sur- 
face layer is based. 



or theoretically via ab-initio simulations. As long as the 
atomistic details affect however only the off-sets of the 
dielectric constant, the electron affinity, and the effective 
mass, the graded interface model already incorporates 
these details by a suitable parameterization. What is not 
well described is the non-universal region a few atomic 
units below and above the crystallographic interface. In 
particular, intrinsic surface states (Shockley and Tamm 
states^) and additional surface states which may arise 
from the short-range surface potential due to impurities, 
imperfections, and ad-layers are not included. If unoccu- 
pied these states could trap the electron surface layer in 
the vicinity of the interface even for surfaces with posi- 
tive electron affinity where image states are absent and 
cannot trap the surplus charge in front of the surface. 

The remaining paper is structured as follows. In Sec.HIl 
we first construct a crude model for the plasma-induced 
electron surface layer at the interface between a plasma 
and a dielectric wall. It does not account for the merging 
of the electron surface layer with the bulk of the dielec- 
tric. As long as the primary interest is in the region close 
to the crystallographic interface and the band gap of the 
dielectric is large enough, the crude model is sufficient. 
Section Hill describes a refinement of the model which en- 
ables one to also investigate the cross-over of the electron 
surface layer to the bulk of the dielectric. This is partic- 
ularly important for dielectrics with small energy gaps. 
Numerical results for the potential and the electron dis- 
tribution are given in Section HVl and a short summary is 
formulated in Sec. IVl 
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II. CRUDE ELECTRON SURFACE LAYER 

As depicted in Fig. [l] we consider an ideal, planar in- 
terface at z = with the dielectric occupying the half- 
space z < and the discharge occupying the half-space 
z > 0. Chemical contamination and structural damage 
due to the burning gas discharge are discarded. At the 
moment we focus on the physical principles controlling 
the electronic properties of the plasma- wall interface. In 
the model we propose the plasma-induced wall charge to 
be treated as an electron surface layer (ESL) which is 
an interface specific electron distribution on top of the 
charge rc-distribution due to the truncation of the solid. 
The ESL is assumed to be thermalized with the solid 
and to stretch from the plasma sheath over the crystal- 
lographic interface to the bulk of the dielectric. 

The boundary between the ESL and the plasma sheath 
is located in front of the surface at z = zq- It is the 
position where the attractive force due to the surface 
potential 4>surf and the repulsive force due to the sheath 
potential (jjsheath balance each other. Thus, zg is given 

by 



(1) 



It gives the position of an effective wall for plasma elec- 
trons and ions at which, for instance, the fiux balance 
condition of the perfect absorber model, je = with je 
and jj, respectively, the electron and ion flux towards the 
dielectric surface, has to be fulfilled. For z < zq an elec- 
tron is attracted to the surface and thus contained in the 
ESL while for z > zq it is repelled back into the plasma. 
On the solid side, for z < 0, the ESL is bounded because 
of the shallow potential well formed by the restoring force 
from the positive charge in the plasma sheath. 

In this section we will outline the essential building 
blocks of the ESL model. Putting together concepts 
from plasma as well as surface physics a detailed, self- 
contained account seems to be helpful. 



A. Plasma sheath 

In the traditional view, electrons missing in the posi- 
tive space charge region in front of the plasma wall accu- 
mulate on the wall and give rise to a wall potential. For 
the construction of our one-dimensional interface model 
we need the total number per unit area of missing sheath 
electrons (that is, the total surface density of missing 
sheath electrons) as a function of the wall potential be- 
cause it is this number of electrons which can be dis- 
tributed across the ESL. Hence, we require a model for 
the plasma sheath. 

For simplicity, we use a coUisionless sheatbi^, more re- 
alistic sheath modela^^rii make no difference in principle. 
In the coUisionless sheath electrons are thermalized, that 
is, the electron density = exp(e(/)//cBTe), with (j) 
the potential, uq the plasma density and Te the electron 



temperature. The ions enter the sheath with a directed 
velocity Vio and satisfy a source-free continuity equation, 
d(niVi)/dz = 0, implying riiVi = rioWjo, and an equation 
of motion M{vi-^Vi) = —e-^(j), with rii the ion density, 
and M the ion mass. The potential cfi satisfies Poisson's 



equation d^^/dz 



~47re(n,; — Uf,). Thus, the governing 



equations for the coUisionless plasma sheath are 



dvi 

dz2 



and 

Mdz 



Wo 

exp 

Vi \kBle 



e(j) 



Using dimensionless variables 



z 



and u 



where 



A 



D 



equations ([2]) and ([3]) become 



and 



M 



Uq 



and 



?; = cxp(-77) 

u 



(2) 
(3) 

(4) 

(5) 



(6) 
(7) 



In the ESL model the plasma occupies not the whole 
half space z > but only the portion z > zq (see 
Fig-IIJ- The integration of the first equation gives 
u = — y^2r] + Uq, where uq — Vio jcg is the reduced ve- 
locity of ions entering the sheath, so that the second 
equation becomes 



- exp(-77) . 



(8) 



Using the boundary condition that the potential and the 
field vanish far inside the plasma, that is, r/ — > and 
if — > for ^ — > oo, Eq. (|5]) can be integrated once and 
we obtain 



-2uQ^2iri + ul + 2 cxp(-7y) + 2uq 



ul-2 



(9) 



For ions entering the sheath with the Bohm velocity 
Uq = —1. The field at the wall as a function of the wall 
potential rj^j ~ ^(^o) is then given by 



= -v/2V27^7 + T+2exp(- 



(10) 



The total surface density of electrons in the ESL equals 
the total surface density of missing sheath electrons, in 
other words, the total surplus surface density of positive 
ions in the sheath N which can be calculated from the 
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electric field at the wall. Integrating Poisson's equation 
yields 



1 d'^0 
dzK-n.) = -— dz — 



1 d0 
47re dz 



(11) 



Combing Eqs. (|lip and (fTO)) gives the total surface den- 
sity of electrons to be inserted into the ESL as a function 
of the wall potential. 

The wall potential itself is determined by the flux bal- 
ance condition, je = ji, which, in the ESL model, is 
assumed to be fulfilled at z = zq. Using the Bohni fiux 
for the ions and the thermal flux for the electrons. 



no 



—— and je = T»^oi 

M 4 V TTTTlp 



the wall potential is given bjM* 



that is, 



1 f M 

2 V ZniTLr- 



In 

2e 



27rmf 



, (12) 



(13) 



(14) 



In the collisionless sheath model the wall potential de- 
pends only on the electron temperature and the ion to 
electron mass ratio. 



B. Surface dipole 

We now turn to the interface region in which the miss- 
ing sheath electrons will be inserted. This region is ab- 
sent in the traditional modeling of plasma walls. In our 
model it is an extended region surrounding an ideal di- 
electric surface. In comparison to the electrons respon- 
sible for the chemical binding within the dielectric the 
additional electrons coming from the plasma are only a 
few. The electronic structure of the surface, in partic- 
ular, the charge re-distribution due to truncation of the 
solid and the offset of the energy bands in the bulk with 
respect to the potential outside the dielectric will not be 
changed significantly by the presence of the surplus elec- 
trons comprising the wall charge. 

In order to quantify the above statement let us first 
consider the electrostatic potential and the electronic 
structure of a free-standing, uncharged dielectric surface. 
According to Tung)^ it has to minimize the thermody- 
namical potential and satisfy Poisson's equation imply- 
ing that the potential is continuous across the surface. 
Strictly speaking, the continuity of the potential only 
applies to the microscopic crystal potential which has 
to merge continuously with the surface potential outside 



the crystal. The averaged long-range potential, in con- 
trast, can be discontinuous at the interface. It is this 
offset which is encoded in the surface dipole. 

The energy of an electronic state in the bulk of the 
dielectric can be referenced to the vacuum level V^(oo) = 
0, that is, the potential far outside the crystal, in the 
following way,— 



Eik{r) = eik - eVceii - eVsir) , 



(15) 



where eik is the quantum-mechanical contribution to the 
energy, Vceii is the averaged potential of a cell due to 
the charge distribution within the same cell, and Vs^r) 
is the long range potential due to the surface dipole, 
space charges, and external fields. In the simple two- 
band model depicted in Fig. [1] i = i;,c. Vsir) contains 
the surface dipole arising from the truncation of the solid 
and responsible for the potential offset at the surface and 
a slowly varying component due to external fields and in- 
ternal and plasma-induced space charges. External fields 
and internal space charges will be neglected in the follow- 
ing and plasma-induced space charges will be accounted 
for by Poisson's equation (see below). 

In order to judge whether the surplus charge arising 
from the plasma affects the surface dipole it is useful to 
consider first the typical strength of the surface dipole of 
a free-standing, uncharged dielectric surface. It results 
from a charge double-layer in immediate proximity to 
the surface. Depending on the material it can have var- 
ious origins. For an ionic crystal, for example, it is the 
lattice relaxation at the surface which makes anions or 
cations to protrude and the other species to retract (e.g. 
protruding oxygen and retracted cations for magnesium 
oxide»2i) , while for semiconductors it is the regrouping of 
covalent bonds which leads to charge re-distribution at 
the surface. Even in the absence of these effects the mini- 
mization of the thermodynamic potential of the surface's 
electrons leads already to an electron density leaking out 
into the vacuum. This is particularly important for met- 
als. As a result a charge double- layer is formed over a 
length on the order of a lattice constant. 

The dipole layer is usually characterized by a dipole 
strength 



eD = eVsir;)-eV{r^) 



(16) 



where Vs {r~ ) is the limit of the long range potential just 
inside the crystal at the surface position and V{f^ ) 
is the limit of the potential at that position just outside 
the surface. Usually these two potentials, which char- 
acterize the discontinuity of the long-range potential at 
the surface, are termed the potential just inside and the 
potential just outside, respectively4^i^ Here, just out- 
side denotes a distance small compared to variations of 
the long-range potential due to external fields or space 
charges but large compared with the width of the charge 
double-layer. Note also that in the definition of the po- 
tential just outside the image potential is assumed to 
have already decreased to zero 4^ This point will be im- 
portant later. 
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The strength of the dipole layer is a microscopic prop- 
erty of the surface which is relatively insensitive against 
the additional charges from the plasma. The reason for 
this lies in the small number of additional electrons from 
the plasma compared to the number of displaced elec- 
trons involved in the formation of the dipole layer. To 
prove this statement wc give a simple estimate. Typical 
surface dipoles eD arc on the order of electron Volts. 
For a double-layer of one A to each side of the crys- 
tallographic interface a potential difference of one Volt 
requires a surface charge density of 5.5 x lO^^cm'^. 
The surface charge density at the wall of a helium dis- 
charge with plasma density hq = lO^cm^"^ and electron 
temperature fc^Tg = 2eV amounts, however, only to 
4.4 X lO^cm"^. For typical plasma densities, the number 
of additional electrons is thus far too small to lead to a 
change of the surface dipole. Even for semiconductor- 
based microdischarges^iiS, which can have much higher 
plasma densities, we expect the surface dipole of the 
plasma wall not to be modified by the plasma. 

In view of the just given estimate, we have to revise an 
assumption in our previous work)2^ where we assumed 
the surface charge accumulating on the wall would in- 
crease the dipole energy eD by e4>w , leading to the image 
states being pushed from the band gap into the energy 
region of the conduction band. The numbers given in 
the previous paragraph indicate, however, that the band 
line-up of the conduction band and the potential just 
outside the solid, will not be affected much by the wall 
charge. Hence, if a negative electron affinity supports im- 
age states in front of the uncharged surface, these states 
remain in the band gap for the charged surface. Electron 
trapping as investigated in Refs i^^i^^ is thus even possible 
for charged plasma walls. 

Instead of the dipole strength eD which cannot be mea- 
sured directly, it is more convenient to characterize the 
dipole layer by the electron affinity x which is a mea- 
surable quantity for a dielectric surface4^ The electron 
affinity is the energy released when an electron is moved 
from just outside the surface to the bottom of the con- 
duction band. It accounts for charge re-distribution in 
the vicinity of the surface due to the truncation of the 
crystal. While many surfaces have positive electron affin- 
ity such as AI2O3 or Si02, there are also materials with 
negative electron affinity, for instance, diamond^, boron 
nitride^, or the alkaline earth oxides<22^ The electron 
affinity depends also on ad-atoms. In some cases this 
is even used to control the electron affinity of a surface. 
Terminating, for instance, a surface with weakly elec- 
tronegative elements such as hydrogen induces a nega- 
tive electron affinity;^ while termination with strongly 
electronegative elements can lead to a positive electron 
aflinity4^ 

From Eq. (jTS]) it is clear that x equals eD plus a bulk 
contribution, 

X = -eVs{f^) - E,k{r7) = eD - Ec + eV,M , (17) 
where Ec — ecm denotes the minimum of the conduc- 



tion band. We can thus use x to characterize the po- 
tential offset at the surface. There is however a caveat. 
The long-range potential inside the solid is only specified 
up to a constant 4^ Typical choices are the cell- averaged 
potential or the inter-sphere potential of the muffin-tin 
approximation. For our purpose it will be however more 
convenient to take the conduction band minimum as the 
long-range potential inside the solid. This choice is mo- 
tivated as follows. We are considering a dielectric with a 
large energy gap. The valence band is thus fully occupied 
and the conduction band is essentially empty. Hence, 
only the conduction band can be populated by additional 
electrons coming from the plasma and referencing the 
electrostatic potential inside the solid to the conduction 
band minimum allows us to relate the total surplus elec- 
tron density in the interface region to the potential in the 
interface region in analogy to what we have done for the 
plasma sheath in the previous subsection. 

Adopting the above discussion to the one-dimensional 
model shown in Fig. [T] and assuming a quadratic disper- 
sion for the conduction band, the energy of an electron 
in the conduction band is given by 

Ekiz) = - ecPsurfiz) , (18) 

where (t>surf{z) is the total surface potential to be deter- 
mined in the next subsection and the offset of the elec- 
trostatic potential at the surface, 

e(f>surf{0^) - e(j>surfiO^) = X , (19) 

encompasses the surface dipole as well as the unspecified 
bulk contribution. 



C. Image potential 

The surface potential of the bare, uncharged surface 
comprises at least the surface dipole and a long-range 
contribution, the image potential, resulting from the mis- 
match of the dielectric constants at the surface. Far away 
from the surface the image potential is given by^ 



But this expression cannot be employed for our purpose 
because the singularity at z = prohibits a smooth elec- 
tron distribution across the interface. In reality the im- 
age potential has to continuously merge with the crystal 
potential. Equation pO[) is thus also unphysical. 

To obtain a realistic image potential without perform- 
ing an atomistically accurate calculation we employ the 
model of a graded interface. It also has the virtue to 
be parameterizable with experimentally measured val- 
ues for the electron affinity, the dielectric mismatch, and 
the mismatch between effective electron masses. The 
model incorporates therefore important properties of a 
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surface, most importantly, it accounts for the charge re- 
distribution due to the truncation of the sohd. 

InitiaUy proposed by Sterol to remove the unphysical 
singularity of the image potential at the interface of two 
dielectrics, the graded interface model assumes the di- 
electric constant e to vary over a distance on the order of 
a lattice constant. Later the model was extended to vari- 
ations of other physical quantities and applied to semi- 
conductor heterostructures and nanopores i^^'^^ Clearly, 
because of the interpolation the model cannot account for 
effects associated with intrinsic surface states (Shockley 
and Tamm states^) and additional surface states which 
may arise from the short-range surface potential. Nev- 
ertheless, the graded interface model is a reasonable de- 
scription of a surface. 

In the spirit of a graded interface, we assume the dielec- 
tric constant e. the electron mass m, and the potential 
offset at the surface to vary smoothly according to the 
grading function 



■ sm I 



z < ~a 
-a < z < a , 
z > a 

(21) 



where a is the half width of the graded interface and 
stands for the quantity that varies across the interface. 
We use the value a = 5 A which is an estimate used in pre- 
vious applications of the graded interfacci^"— While the 
grading parameter a is not based on definite experimental 
or theoretical results it is motivated by the assumption 
that the bonding and the electron density at the surface 
change over one to two lattice constants implying a tran- 
sition layer for the effective potential and the dielectric 
constant that is somewhat larger. Hence, across the in- 
terface the electron mass, the dielectric constant and the 
offset potential are given by 



fn{z) = gm.},,m,Az) , <^{z) = 9es{^) ' (22) 



and 



<t>offset{z) = -9xfl{z) 



(23) 



respectively, with the effective mass of the conduc- 
tion band. 

Within the model of the graded interface, the image 
potential is the change in the selfcnergy of an electron due 
to the proximity of the dielectric mismatch. Positioning 
the electron at tq it is given bj*^ 

0.™(n)) = ^[0"(r'o)-0°(ro)] , (24) 

where 0™ (f) is the potential in the medium with dielec- 
tric mismatch arising from the electron at ro and (ffiif) 
is the same quantity in a homogeneous medium with di- 
electric constant e(zo). Hence, </)"'(r) is the solution of 



V (e(0)V(?!."(f)) = 47re,5 (f - ro) 



(25) 



while (tP{r) is the solution of 

47re 



V^0"(f) 



e(^o) 



&{r- ro) 



(26) 



To solve Eqs. ([25]) and (HH) we follow Stern^^ and make 
the ansatz 



</>™'0(z,p. 



1 °° 

— y 



dqqJi{pq)Ji{poq) 



X e''('^-'^'''A™'°(z) . (27) 



Placing the electron on the z axis, po ~ which implies 
Jq{Pqq) = 1 E^nd Ji{poq) = for I > 0. Hence, we need 
to keep only the I ~ term, so that 

1 r°° 

cj,"^^\z,p,4,)= dggJo(p<z)A;"'0(z) , (28) 

where A^^{z) is given by 



-27re 



-q\z-Zo\ 



(29) 



e(2:o)g 

and A'^{z) is the solution of 

A™"(z) + ^A™'(z) - q^A^{z) = ^Siz - zo) 
^ e[z) ^ ^ e(z) 

(30) 

which has to be obtained numerically. The image poten- 
tial is then given by 



<t'im.{zo) 



1 

47r 



dOT(^;"(zo)"AO(zo)) . (31) 



In contrast to Eq. (|20)) it is now smoothly varying across 
the interface with a deep well on the low-e side and a 
small bump on the high-e side. 

The total surface potential comprises the graded offset 
potential p3)) and the graded image potential. Hence, 



(t'surfiz) = 4>im{z) + 4>offset{z) 



(32) 



It is continuous across the crystallographic interface at 
z = and enables us thereby to also calculate a smoothly 
varying electron distribution in the ESL. The band struc- 
ture and the total surface potential at the graded inter- 
face arc visualized in the lower panel of Fig. [1] 

Using Eq. ([1]) we can now determine the position zo of 
the effective wall, that is, the maximum extent of the ESL 
on the plasma side. The derivative of the bare surface 
potential is (t>'aurf — ^'offset + ^'im- Due to the relatively 
weak field in the sheath compared to the image force, 
the boundary zq will be so far away from the interface 
that 4>'offset vanishes and the image potential obeys ^U^. 
Thus, the boundary between the ESL and the plasma 
sheath is given by 



Zo 



{e-l)e 



, 4(e+l)0^ 

with 0'^ = ^{kBTev'J/ieXD) and t]'^ given by 1^. 



(33) 
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D. Electron distribution 

The plasma-induced wall charge is assumed to be in 
thermal equilibrium with the wall. Hence, the distribu- 
tion of the excess electrons in the ESL has to minimize 
the excess electron's grand canonical potential in the ex- 
ternal potential due to the surface. The coupling to the 
sheath is maintained by the constraint that only as many 
electrons can be filled into the ESL as are missing in the 
sheath and the boundary conditions to the Poisson equa- 
tion which links the electron distribution in the ESL to 
the (internal) electrostatic potential. 

To minimize the grand canonical potential of the sur- 
plus electrons wc follow Tkharev and Danilyuk^ and 
apply density functional theory^ii^ to the graded inter- 
face. While more refined schemes of density functional 
theory^ could, in principle, be employed, we will use for 
the purpose of this exploratory calculation density func- 
tional theory in the local approximation. Quite generally, 
the grand canonical potential of an electron system in an 
external potential V{'r) is given in the local approxima- 
tion by 



n — I V{r)n{r)dr— — j 0c'(r)n(r)df 



+ G[n] — H n{r)df 



(34) 



where G[n\ is the grand canonical potential of the ho- 
mogeneous system with density n{r) and the Coulomb 
potential is determined by 



V (e(f)V0c(^)) = 47re?i(r). 



(35) 



The ground state electron density minimizes fl, that is, 
it satisfies 



V{f} - e(t)c{r) + y!'{n) - ^ = 



(36) 



where ijp'{n) — 5G[n]/5n is the chemical potential for the 
homogeneous system. 

Specifically for the excess electrons in the one- 
dimensional graded interface Eq. (|36p reduces to 



'e4){z) -f ^''(z) - /i = 



(37) 



where ij.^{z) = fi''-{n{z),T) and the electrostatic poten- 
tial. 



(2) = 4>surfiz) + 4>ciz) 



(38) 



consists of the potential of the bare surface given by 
Eq. (fT9|) and the internal Coulomb potential which sat- 
isfies Poisson's equation. 



'^iz)-T-4'ciz) ) = 47ren(z) , 



(39) 



with the graded dielectric constant e{z) given by Eq 
and the boundary conditions (j)c{zo) = (j)w and 0p(zo) 



(j)'^ to guarantee continuity of the potential at and to 
include the restoring force from the positive charge in the 
sheath. Note that the Coulomb potential derived from 
this equation includes the attraction of an electron to the 
image of the charge distribution. 

For the functional relation /i'*(z) = fi^{Ti{z),T) we 
take the expression adequate for a homogeneous, non- 
interacting, non-degenerate electron gas. 



n{z) 



1 / m{z)kT 

71 



(40) 



This is justified because the density of the excess elec- 
trons is rather low and the temperature of the surface is 
rather high, typically a few hundred Kelvins. 

In order to calculate the quasi-stationary distribu- 
tion of the surplus electrons, Eqs. (|37|) and (p9)) have 
to be solved sclf-consistently with the additional con- 
straint that the total electron surface density in the ESL 
equals the total surface density of electrons missing in 
the sheath, that is. 



dzn{z) = TV 



(41) 



with N given by Eq. (|lip . In the above equation we in- 
troduced a cut-off Zs < at which the ESL terminates in- 
side the dielectric. As long as \zs\ is chosen large enough 
it does not affect the numerical results close to the sur- 
face. An improved treatment of the ESL, avoiding the 
ad-hoc cut-off, is given in the next section. 

Within the crude ESL model developed in this section 
the computation is performed iteratively in the interval 
Zs < z < zq according to the following scheme: (i) We 
start with the potential of the empty surface given by 
Eq. (inSl) with 4)c{z) obtained from Eq. with n{z) = 
but with the boundary conditions at zq as specified, (ii) 
We integrate both sides of Eq. (PD|) over z with fi^iz) 
given by Eq. (jST)) . Enforcing the constraint (|4T|) deter- 
mines /i. (iii) Using /x we calculate from Eq. (|37p a new 
fj.^{z) which gives with Eq. (|^(])) a new electron density 
n{z). (iv) Lastly, we determine from Eq. (|39p the electro- 
static potential associated with the up-dated n{z). Steps 
(ii) - (iv) arc iterated until /i, which is far below the con- 
duction band edge because of the non-degeneracy of the 
excess electrons, converges. 



III. REFINED ELECTRON SURFACE LAYER 

In the previous section wc have taken into account only 
the electron concentration in the conduction band of the 
dielectric due to the electrons coming from the plasma. 
For wide band gap materials this is justified, especially 
near the surface, as their concentration is much larger 
than the intrinsic carrier concentration. Deep inside the 
dielectric, however, charge neutrality is not enabled by 
a vanishing electron density, but by the electron density 
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FIG. 2: Sketch of the refined model of the interface between 
a plasma and a dielectric wall. In the plasma equal densities 
of electrons and ions ensure quasi-neutrality. The positive 
space charge in front of the effective wall defines the plasma 
sheath. The ESL contains a very narrow interface-specific re- 
gion (ISR) where the model of the graded interface is used 
and a wide space charge region (SCR) which allows a con- 
tinuous merging with the neutral bulk of the dielectric where 
intrinsic electrons and holes balance each other to guarantee 
charge neutrality. Note, the widths of the various regions are 
not to scale. 



decreasing to its intrinsic value, which is then balanced 
by the intrinsic hole concentration in the valence band. 

To take this effect into account, vifhich is particularly 
important vifhen the additional electrons accumulate deep 
inside the bulk of the dielectric, we divide the ESL into 
two regions: a very narrow interface-specific region (ISR) 
and a wide space charge region (SCR) in the bulk of the 
dielectric. The parameter denotes now no longer an 
ad-hoc cut-off but the boundary between the two regions. 
It has to be chosen so that the ISR includes the major 
effect of the image potential in the dielectric implying 
Zs < —zq. The electron distribution and the potential in 
the ISR are calculated using the density functional ap- 
proach outlined in the previous section. In the SCR we 
use for simplicity the model of an intrinsic semiconduc- 
tor to describe electron and hole densities as well as the 
long-range potential. As the energy bands in the dielec- 
tric follow the long-range potential the refined ESL also 
captures the band bending which might be induced by 
the presence of the wall charge. It is however only sig- 
nificant when most of the excess electrons are trapped in 
the SCR and not in the ISR. 

Figure[5]schematically shows the electron and hole den- 
sities for the refined ESL model. The boundary between 
the plasma sheath and the ESL is still located at zq. As 
our model does not encompass the electron and ion flux 
from the plasma for z < zq, the densities rig and Ui are 
discontinuous at zq. This is obvious for the ions which 
are not allowed to enter the solid. The discontinuity of 
the electron density, in contrast, arises because we en- 
sure only the total number of missing sheath electrons 
per unit area to be conserved. This global constraint 
cannot guarantee continuity of the electron density at 
Zq. At the boundary between the ISR (zs < z < zq) and 
the SCR {z < Zs) the electron density and the poten- 



tial are continuous. In principle, also the hole density p 
should be continuous. As p{zs) <C n(zs) for the materials 
we are considering, we can however neglect holes in the 
ISR. 

For the modeling of the SCR it is convenient to use 
ipiz) = (f>{z) — 4'huik for the long range potential, which 
vanishes for charge neutrality in the bulk. Here, ^buik = 
0(— oo) (see below for an explicit relation for 4>huik)- 
Then, Poisson's equation is given by. 



dSp{z) 

dz2 



An 



(— en(z) + ep{z)) 



(42) 



where the electron and hole densities for an intrinsic semi- 
conductor with parabolic bands whose extremal points 
are, respectively, Ec and Ey are given bj*^ 



p{z) = 




V2\ TTh^ 



(43) 



(44) 



From a comparison of the exponents in Eq. (j43p and 
Eq. (PO]) , where fi^ is given by Eq. ([57)) we find 



jy = ^1 + Ec + e(t>buik 



(45) 



Far from the surface, ip — Q and n = p = nt. This gives 
the chemical potential 



Ey+Ec , 3,^ ^, (m^ 
V = V -kBTlvL 



and the bulk carrier concentration rif, 



rifa 



(46) 



(47) 



where E'g = Eq—Ev- Hence, Poisson's equation becomes 



d'^^p{z) Aire / £^ 



d22 



47re / <"»'^) 



and using dimensionless variables 

eip , . z- Zs 

V = 1 — 7^ sind t, 



knT 



Ld 



with Ljj = ^ eksT / Ane'^ni, we obtain 



7j" = e" - e^" 



(48) 



(49) 



(50) 



This equation can be integrated once, which gives 

(7/)2 =4cosh(77)-FC . (51) 

The boundary conditions in the bulk rj ^ and rj' ~ 
for ^ — oo imply C = —4 so that Eq. ([ST]) becomes 



77' = V8sinh(| 



(52) 
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Integration with the boundary condition at Zg, that is, at 
^ = 0, 77(0) = 77^ and requiring 77 — > for ^ -> —00 gives 



77±(0=T21n 



± tanh 



with 



±2artanh 



exp 



V2 



TVs 
2 



(53) 



(54) 



where the upper sign is for rjs > and the lower sign for 
Vs < 0. 

In analogy to what we have done at the boundary of 
the ESL with the plasma sheath at z = zq we relate the 
potential rjs to the total electron surface density in the 
space charge region. From Poisson's equation we obtain 
for the total electron surface density in the SCR 



TV' 



SCR 



2 „ d?7 



dz{n — p) = LjjUf, 



dz 



where 77' is given by Eq. ([5^ . so that 

/ jySCR 



(55) 
(56) 



or. 



rjs — 2arsinh 



\V8L 



Drib 



(57) 



For a negative space charge rjs > 0, so that the poten- 
tial is given by ip{z) = {ksT / ej-q'^ {{z — Zs)/Ld) and 
the electron and hole densities are given by n{z) = 
rn^(fi+((^-'^s)/Ln) and p{z) = nbe-'^^^^^-^^^/^"\ respec- 
tively. The relation between ip and (f> is given by 
(kbuik = 0(z) - ipiz). Since ip{zs) = {kT/e)ris we obtain 
(kbuik = (pizs) - {kT/e)r]s. 

Now, quite generally, the excess electrons in the ESL 
are distributed over the ISR and SCR according to 



SCR/ 



(58) 



where N^^^ is the surface density of electrons in the ISR, 
H is the chemical potential in both regions, and N is the 
total surface density of missing sheath electrons given by 
Eq. pT|) . The total surface density in the ISR is given 
by N^^'^ ~ J^" dzn{z), where n{z) is calculated with 
the density functional approach for the graded interface. 
Requiring continuity of the electron density at z^, 



-(- 



■Kh? ) 



e ''B 



,(/i+e0(z,)) 



me ' 



(59) 



gives 77s as a function of fi. From 77^ we finally obtain 
using Eq. §6^ N'^^^ifj,). 

For the calculation of the electron distribution and the 
potential in the refined ESL we use the iteration cycle de- 
scribed in the last section with one modification. In step 



TABLE I: Material parameters for the dielectrics considered 
in this work: dielectric constant e^, electron affinity Xi con- 
duction band effective mass nic, valence band effective mass 
771^, band gap Eg. 

X [eV] 771g [771e] my [me] Eg [eV] 

MgO 9.8 [51J -0.4 \28} 0.4 [52] 

AI2O3 9.9 [53| 2.5 [541 0.4 [55] 4.0 [55] 8.8 [56] 

Si02 3.78 [57] 1.3 [54] 0.5 [58| 0.58 [59] 9.2 [541 

GaAs 13.1 [60] 4.07 [60] 0.067 [60] 0.45 [60] 1.42 [60| 



(ii) we solve Eq. instead of Eq. (pij) to fix /i. From 
fi we obtain using Eq. (|59p 7]^ which in turn determines 
the electron distribution and the potential in the SCR. 
This gives for each iteration step a continuous potential 
and electron distribution at z^. As before, the steps (ii) 
-(iv) are iterated until fj, converges. 

At the end of this section let us finally mention two 
simplifications of the refined ESL model which could be 
used, respectively, for large band gap dielectrics irrespec- 
tive of the electron affinity and dielectrics with small 
band gap and positive electron affinity. In the former 
case the intrinsic carrier concentration Ub is very low and 
the merger with the bulk occurs very deep in the di- 
electric. Almost all surplus electrons are however much 
closer to the surface where the holes can be neglected. 
This can be seen from the differential equation for 77. For 
small Ub Eq. ((57|) gives a large rjs. As 77 satisfies a highly 
nonlinear differential equation (|50p . a large 77s implies a 
steeper gradient of 77 near the surface so that almost all 
electrons are concentrated close to the surface where ne- 
glecting the holes has little effect. Hence, for large band 
gap dielectrics surplus electrons can be filled into a suf- 
ficiently large ISR for which the crude ESL model of the 
previous section will be sufficient provided the cut-off Zs 
is large enough. The merger with the bulk is of course 
not correctly captured by such an approach. 

For dielectrics or semiconductors with small energy 
gaps and positive electron affinity, on the other hand, 
almost all surplus electrons are deep inside the material. 
It is thus a good approximation to neglect the ISR and to 
fill all electrons in a SCR. Neglecting the surface potential 
has little effect in this case and using the SCR already for 
z <0 gives a good description of the electron distribution 
inside the ESL. The electron density and potential at the 
surface and in front of it can of course not be captured 
by such an approach. As before ip{z) = (j){z) — (pbuik with 
4>buik = 4>{zs) ~ {ksT / e)rjs where (l){zs) is now the limit 
of the long range potential just inside the dielectric given 
by (j){zs) = +x/e- 



IV. RESULTS 

We now use the ESL model to calculate for a helium 
discharge in contact with a MgO, AI2O3, and Si02 sur- 
face the potential and the density of excess electrons 
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FIG. 3: Plasma-supplied excess electron density n (upper 
panel) and the potential — ^ (lower panel) it gives rise to 
for a MgO surface in contact with a helium discharge with 
no — W^cm^'^ and fcsTe = 2eV calculated without account- 
ing for a SCR in the dielectric (crude ESL model). The cut-off 
of the interface region is Zs = —zo- As can be seen, almost 
all of the plasma-induced wall charge is located in the well of 
the image potential in front of the surface. 



across the plasma wall. Our main focus lies in the 
identification of generic types of electron distributions 
in the ESL depending on plasma and surface param- 
eters. Unless otherwise stated, we use for the plasma 
density uq = lO^cm""^ and for the electron temperature 
ksTe = 2eV. The parameters of the dielectric surfaces 
are given in Table ID Preferentially we used experimental 
data for the various quantities only if not available wc 
employed theoretical valuesi ^^'^^" — 

First we give typical values for zq, the position where 
the ESL merges with the plasma sheath. It is calculated 
from Eq. ([T|) and should thus depend not only on plasma 
but also on surface parameters. Our results for MgO 
{zq 6.08 X 10-''^cm), AI3O3 (zq = 6.09 x lO'^cm), 
and Si02 (20 = 5.14 x 10~^cto) indicate however that 
Zq is relatively insensitive to e which is the only surface 
parameter affecting zq when the sheath is assumed to 
be coUisionless. Even the significantly smaller e of Si02 
docs not alter zq considerably. For the helium discharge 
considered zq is irrespective of the dielectric always on 
the order of a micron. 

Of particular importance for the distribution of the 
excess electrons in the ESL is the electron affinity x, 
characterizing the offset of the conduction band to the 
potential just outside. For x < (MgO) the conduction 
band minimum lies above the potential just outside. It 
is thus energetically favorable for electrons to be located 
in the image potential in front of the surface. Figure [3] 
showing the electron density and the potential in the ESL 
of MgO verifies this. The energy of an electron in the im- 
age potential —ecj) indeed reaches a minimum just in front 
of the surface at the beginning of the graded interface. 
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FIG. 4: Plasma-supplied excess electron density n (upper 
panel) and the potential —0 (lower panel) it gives rise to for 
an AI2O3 surface in contact with a helium discharge with 
no — 10^cm~^ and ksTe — 2eV. The red lines show data ob- 
tained from the refined ESL model accounting for an ISR and 
a SCR, the boundary between the two was put at Zs ~ —3^0, 
the green circles show data for a model which consists only 
of an ISR with cut-off Zs — —0.9 cm (crude ESL model), 
and the blue triangles show data for a model consisting only 
of a SCR for z < 0. Irrespective of the approximation, the 
plasma-induced wall charge extends deep into the bulk. Note 
the different scales of the axes for the left and right panels. 
On the spatial scale of the SCR shown in the left panels the 
ISR of the right panels becomes a vertical line. 



For negative electron affinity, the excess electrons coming 
from the plasma thus form an external surface charge in 
the image potential in front of the crystallographic inter- 
face. The band bending associated with it is negligible. 
The external surface charge is very narrow, it can thus be 
considered as a quasi-two-dimensional electron gas, sim- 
ilar to the surface plasma anticipated by Emeleus and 
Coulter J^i^ 

For X > 0, on the other hand, the conduction band 
minimum is below the potential just outside. It is thus 
energetically favorable for electrons to accumulate inside 
the dielectric. This can be seen in Fig. U which shows the 
electron density and the potential in the refined (red line) 
and simplified ESL (open green circles and blue triangles) 
for an AI2O3 surface. The surface potential consists of 
an attractive well in front of the surface but the mini- 
mum potential energy for electrons —ecj) is reached inside 
the dielectric. Excess electrons coming from the plasma 
are thus mostly located inside the wall and the electron 
distribution extends deep into the bulk. This extended 
negative space charge also leads to a band bending near 
the surface. Note the different scales of the axes for the 
left and right panels of Fig. 2] On the scale where vari- 
ations in the SCR are noticeable the ISR is basically a 
vertical line. 

If one neglects the SCR and fills all excess electrons 
into the ISR (crude ESL) the potential and the electron 
distribution are correctly described at and close to the 
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FIG. 5: Electron density n and hole density p at a GaAs sur- 
face in contact with a helium discharge with no — VS' cmT'^ 
and fcaTe = 2el/ calculated with the refined ESL model with- 
out ISR. The plasma-induced wall charge sits inside the GaAs 
wall. Deep inside the bulk charge neutrality is achieved by an 
equal density of electrons and holes. 
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FIG. 6: Center of gravity z of the plasma-supplied excess 
electron distribution at a MgO surface (upper panel) and the 
290% value for the electron distribution at an AI2O3 and a 
Si02 surface (lower panel), all in contact with a helium dis- 
charge with no = lO^cm"^ and fcsTe = 2eV ^ as a function 
of the surface temperature Ts- The data shown in the upper 
and lower panel were obtained, respectively, from the crude 
ESL model and the refined ESL model without an ISR. 



surface but not far inside the dielectric (open green cir- 
cles) because the ad-hoc cut-ofF Zg of the crude ESL leads 
to an unphysical pile- up of electrons near Zs . Hence, only 
if Zs is large enough does the crude ESL model give reli- 
able results for the electron density and potential in the 
vicinity of the surface. Filling all electrons in the SCR, on 
the other hand, cannot describe the immediate vicinity 
of the surface correctly which is however on the scale of 
the SCR an infinitesimally narrow region. It gives only 
for z < ~zq a good description, that is, in the region 
where for x > indeed most of the electrons are located 
(blue triangles). 

While the crude ESL model containing only an ISR 
gives the correct electron density near the surface pro- 
vided Zs is large enough, the merger of the ESL with the 
bulk can only be described with the refined ESL model 
including the SCR. This is particularly relevant for ma- 
terials with smaller band gaps and larger intrinsic carrier 
concentrations than MgO, AI2O3, and Si02. To exem- 
plify this we show in Fig.[5]the electron and hole densities 
(upper panel) as well as the potential —'0 (lower panel) 
for a GaAs plasma wall, calculated for an ESL containing 
only a SCR. At the surface the electron density is about 
three orders in magnitude larger than the hole concentra- 
tion. Thus, the gas phase plasma offers the possibility to 
manipulate the electron-hole plasmas by controlling the 
charge carrier density - tantamount to doping - in the 
near surface region of a semiconductor. Deep inside the 
material electron and hole concentrations are equal, lead- 
ing to charge neutrality and a constant potential. The 
band bending due to the extended space charge in the 
ESL is about 0.09 eV. 

Our results for the electron and hole densities and the 
potential in the dielectric depend of course on the model 



for the SCR. We have used for simplicity the model of 
an intrinsic semiconductor which is appropriate for an 
undoped semiconductor without impurities. Depending 
on doping or impurities a variety of models^ could be 
used to take material specific aspects into account. In 
our exploratory calculation we obtain a rather wide SCR. 
Including the effect of impurities, acting as trapping sites 
in the band gap, would probably reduce the depth of the 
SCR considerably. 

To summarize our results up to this point, we find 
that for negative electron affinity the plasma-induced 
electronic surface charge is located in front of the sur- 
face forming a quasi-two-dimensional electron gas, while 
for positive electron affinity the surplus electrons form 
a space-charge layer in the dielectric leading to a small 
bending of the energy bands. 

The two distinct types of charge distributions in the 
ESL are also reflected in the dependence of the width of 
the plasma-supplied electron distribution on the surface 
temperature. Figure [6] shows the center of gravity z of 
the electron distribution for the MgO surface {x < 0) 
and the 2:90% value for the surfaces of AI2O3 and Si02 
(x > 0), where zqq% is implicitly defined by 

/ dz{n{z) -p{z)) = 0.9 N . (60) 

We use the 299% value because it captures the depth of 
the SCR better than z which depends too strongly on 
the few electrons that penetrate very deep into the bulk. 

For negative electron affinity (MgO, shown in the up- 
per panel of Fig.|6]) the external surface charge is strongly 
trapped in the deep image potential so that z changes 
very little with surface temperature. The width of the 
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internal surface charge for dielectrics with positive elec- 
tron afhnity (AI2O3 and Si02, lower panel), however, in- 
creases dramatically with surface temperature. This can 
be understood as follows. The restoring force from the 
positive ions in the sheath binds internal surface charges 
only weakly to the surface. With increasing surface tem- 
perature, however, high-lying states in the conduction 
band get more and more populated. Hence, some elec- 
trons have rather high kinetic energies, are thus less con- 
fined near the surface by the weak restoring force, and 
penetrate therefore deeper into the bulk. As a result, the 
2:90% value decreases strongly with surface temperature. 

Let us now turn to the discussion of the influence of 
the electron temperature kgT^ and the plasma density 
no on the properties of the ESL. These two parameters 
enter through the total surface density of electrons N 
depleting the sheath and accumulating in the ESL. How 
ksTe and uq affect the interface depends therefore on 
the sheath model and the model used for the interaction 
between plasma particles and the surface. For simplicity 
we have used a collisionless sheath model and assumed 
the surface to be a perfect absorber for plasma electrons 
and ions. The results for the properties of the ESL as a 
function of the plasma parameters are thus to be taken 
only indicative. 

The effect of a variation of tiq and ksT^ is most signif- 
icant for surfaces with positive electron affinity. Table [TT] 
shows the effect of the plasma density ng for a Si02 sur- 
face. If no increases, the boundary zq between sheath and 
ESL moves closer to the surface. This, however, does not 
affect the charge distribution much as most of the elec- 
trons occupy the SCR inside the dielectric (as shown in 
Fig- m for AI2O3). More important is that an increase in 
no leads to an increase of the total surface electron den- 
sity N. This entails a stronger restoring force from the 
plasma sheath so that the potential well confining the 
space charge inside the dielectric becomes steeper and 
the electrons in the SCR of the ESL are shifted towards 
the surface, in other words, the 2:90% value increases with 
uq. Mathematically, the steeple-like shape of the elec- 
tron distribution arises because a larger N leads through 
Eq. ((57|) to a larger rjs which makes the potential steeper 
at the surface so that the electron distribution is more 
peaked there. This trend can be seen in the lower panel 
for Fig. H 
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FIG. 7: Plasma-supplied surplus electron density n at a Si02 
surface in contact with a helium discharge as a function of the 
electron temperature (upper panel, no = lO'^ cm~^) and the 
plasma density (lower panel, fcsTe = 2eV) for Ts = 300K. 
The refined ESL model without an ISR was employed to pro- 
duce the data. 



A variation of the electron temperature ksT^ has simi- 
lar effects as the variation of the plasma density. If ksT^ 
increases, the total surface density of electrons increases 
also, as can be seen from Table IIIII As shown in the 
upper panel of Fig. [71 this leads again to a steeple-like 
electron distribution which is the more concentrated at 
the surface the higher the electron temperature is. 

For a surface with negative electron affinity (MgO) the 
surplus electrons are strongly bound in the image poten- 
tial. While a variation of /csTe or no changes the total 
number of surplus electrons per unit area in the same 
way as for a surface with positive electron affinity, the 
distribution of the surplus electrons within the ESL is 
not affected significantly because of the strong image in- 
teraction. 

So far we have shown the potential and the electron 
distribution in the ESL. Now, we will compare potential 
and charge distribution in the ESL with the ones in the 



TABLE II: Surface density of electrons in the ESL A*', wall 
potential (j>w, plasma sheath - ESL boundary zq and the 290% 
value for Si02 in contact with a helium discharge with ksTe — 
2eV for different values of the plasma density no . 



TABLE III: Surface density of electrons in the ESL A'', wall 
potential (j>w , plasma sheath - ESL boundary zq and the 290% 
value for Si02 in contact with a helium discharge at no = 
W^cm~^ for different values of the electron temperature ksTe- 



no [10^ cm~^] 


A' [10^ cm~2] 


<^iv [V] 


zo [lO"'^ cm] 


^90% [cm] 


ksTe 


[eV] N [10^ cm~^] 


<t}w [V] 


Zo [10 ^cm] 


Z(jo% [cm] 


10 


4.38 


-7.07 


5.14 


-0.222 


0.5 


2.19 


-1.77 


7.27 


-0.444 


20 


6.20 


-7.07 


4.32 


-0.157 


1 


3.10 


-3.53 


6.11 


-0.314 


50 


9.80 


-7.07 


3.44 


-0.099 


2 


4.38 


-7.07 


5.14 


-0.222 


100 


13.9 


-7.07 


2.89 


-0.070 


5 


6.93 


-17.7 


4.086 


-0.140 
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FIG. 8: Density of plasma-supplied surplus electrons trapped 
in the ESL as well as electron and ion density in the plasma 
sheath (upper panel) and potential (lower panel) for a MgO 
surface in contact with a helium discharge (no = lO^cm"^ 
and fcsTe = 2eV). The data were obtained from the crude 
ESL model. 
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FIG. 9: Density of plasma-supplied excess electrons in the 
ESL as well as electron and ion density in the plasma sheath 
(upper panel) and potential (lower panel) for a Si02 surface in 
contact with a helium discharge (no = lO^cm^'^ and ksTe = 
2eV). The refined ESL model without an ISR was employed 
to produce the data. 



plasma sheath. The electron distribution at the interface 
is the quasi-stationary electron gas on top of the charge 
re-distribution due to the truncation of the solid that 
guarantees flux equality at the sheath-ESL boundary zq. 
As already mentioned not included in this simple model 
is the flux of plasma electrons and ions in the ESL before 
the electrons are trapped at the surface and the ions rc- 
combine with the negative wall charge. The electron and 
ion densities in this model arc thus discontinuous at zq. 
The potential, however, which has been obtained from 
the integration of Poisson's equation is continuous and 
differentiable everywhere. Between the well of the image 
potential and zq the electron and ion flux from the sheath 
would be important. The neglect of the charge densities 
associated with these fluxes does however not affect the 
potential because they are too small to have a significant 
effect. 

Figure [8] shows the ESL and the plasma sheath in front 
of a MgO surface. Due to the negative electron affinity, 
the plasma-supplied surface electrons are bound by the 
image potential in front of the surface. In Fig. |S1 we 
plot the electron and ion density (upper panel), as well 
as the electric potential (lower panel) over the distance 
from the surface z. Far from the surface, the potential 
approaches the bulk plasma value chosen to be zero. In 
the sheath the potential develops a Coulomb barrier and 
reaches the wall potential at zq, the distance where 
the sheath merges with the surface layer (vertical dot- 
ted line). The wall potential is the potential just outside 
to which the energies of the bulk states are referenced. 
Closer to the surface the potential follows the attractive 
image potential while at the surface the repulsive po- 
tential due to the negative electron affinity prevents the 



electron from entering the dielectric (only scarcely seen 
on the scale of the figure) . 

In Fig. [S] we finally plot the electron and ion densi- 
ties (upper panel) as well as the electric potential (lower 
panel) for Si02. Note the linear z-axis in contrast to the 
logarithmic z-axis of Fig. |S1 Due to the positive electron 
affinity, the excess electrons constituting the wall charge 
penetrate deep into the dielectric and occupy therefore 
the SCR of the ESL. Compared to the variation of the 
electric potential in the sheath the band bending in the 
dielectric induced by the wall charge is rather small as 
indicated by the variation of (j) inside the dielectric. This 
is because e is large and the width of the SCR is narrow 
on the scale of the sheath. Only on the scale of the ISR 
(a vertical line at z = 0) the SCR of the ESL is wide. 



V. CONCLUSIONS 

We have studied the potential and the charge distri- 
bution across the interface of a plasma and a dielectric 
wall treating the plasma-induced wall charge as a quasi- 
stationary electron gas trapped by and in thermal equi- 
librium with the dielectric. Our approach is based on a 
model for a graded surface including the offset between 
the potential just outside the dielectric and the conduc- 
tion band minimum arising from the re-distribution of 
charge due to the truncation of the solid as well as the 
image potential due to the dielectric mismatch at the 
boundary. The missing electrons from the sheath popu- 
late the interface potential and thereby form an electron 
surface layer (ESL) which minimizes the grand canonical 
potential of wall-thermalized excess electrons and satis- 
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fics Poisson's equation. 

Within this model the boundary between the plasma 
sheath and the ESL is given by the distance from the 
crystallographic interface where the potential for the ex- 
cess electrons turns from the repulsive sheath potential 
into the attractive surface potential. This distance is typ- 
ically on the order of a micron. It gives the position of 
an effective wall for plasma electrons and ions and thus 
the portion of the ESL which lays in front of the surface. 
Most of the surplus electrons trapped in the ESL, that is, 
the plasma-induced wall charge, will be, however, much 
closer to the surface or even inside the dielectric depend- 
ing on the electron affinity. 

We presented numerical results for the potential and 
the distribution of the plasma-supplied surplus electrons 
at the interface between a helium discharge and the sur- 
faces of MgO, AI2O3, and Si02, respectively. The elec- 
tron distribution within the ESL strongly depends on the 
electron affinity. For negative electron afBnity, the con- 
duction band minimum is above the potential just out- 
side the dielectric. Hence, it is energetically unfavorable 
for electrons to penetrate into the bulk and the surface 
electrons are bound in the image potential in front of the 
surface. In this case, their spatial profile changes little 
over a variation of the surface temperature or the plasma 
parameters. For positive electron affinity the conduction 
band minimum is below the potential just outside the di- 
electric and the surface-bound electrons accumulate in- 
side the wall. The space charge in the bulk broadens if 
the surface temperature is increased and becomes more 
peaked if the total surface density of the electrons miss- 
ing in the sheath is raised through either an increase in 
the plasma density or the electron temperature. 



Separating the ESL into an interface specific and a 
space charge region and modeling the bulk of the dielec- 
tric as an intrinsic semiconductor we also investigated 
how the ESL merges with the bulk of the dielectric. This 
is particularly important for dielectrics with small en- 
ergy gaps and positive electron affinities where excess 
electrons coming from the sheath accumulate not in the 
image potential in front of the surface but deep inside 
the wall. In this case the wall charge may also induce a 
significant band bending. 

Whereas the crude ESL model we proposed neglects 
the space charge deep inside the bulk of the wall and 
is thus only applicable to large band gap dielectrics 
with negative electron affinity where basically the whole 
plasma-induced wall charge is trapped in the image po- 
tential in front of the surface, the refined ESL model 
keeping the interface specific as well as the space charge 
region of the ESL provides a quantitative description 
of the whole spatial structure of the extended charge 
double-layer which forms at a dielectric plasma wall as a 
result of the electrons in the ESL and the positive space 
charge in the plasma sheath. 

The ESL can be regarded as that part of the plasma 
sheath which is inside the plasma wall. It is thus the 
ultimate boundary of a bounded gas discharge and con- 
stitutes, depending on the electron affinity either a quasi 
two-dimensional electron plasma in front of the wall or 
an electron(-hole) plasma inside the wall. 
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